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Abstract 

Background: Whether myofibers increase with a pulsed-wave mode at particular developmental stages or whether 
they augment evenly across developmental stages in large mammals is unclear. Additionally, the molecular 
mechanisms of myostatin in myofiber hyperplasia at the fetal stage in sheep remain unknown. Using the first 
specialized transcriptome-wide sheep oligo DNA microarray and histological methods, we investigated the gene 
expression profile and histological characteristics of developing fetal ovine longissimus muscle in Texel sheep (high 
muscle and low fat), as a myostatin model of natural mutation, and Ujumqin sheep (low muscle and high fat). 
Fetal skeletal muscles were sampled at 70, 85, 100, 120, and 135 d of gestation. 

Results: Myofiber number increased sharply with a pulsed-wave mode at certain developmental stages but was 
not augmented evenly across developmental stages in fetal sheep. The surges in myofiber hyperplasia occurred at 
85 and 120 d in Texel sheep, whereas a unique proliferative surge appeared at 100 d in Ujumqin sheep. Analysis of 
the microarray demonstrated that immune and hematological systems' development and function, lipid 
metabolism, and cell communication were the biological functions that were most differentially expressed 
between Texel and Ujumqin sheep during muscle development. Pathways associated with myogenesis and the 
proliferation of myoblasts, such as calcium signaling, chemokine (C-X-C motif) receptor 4 signaling, and vascular 
endothelial growth factor signaling, were affected significantly at specific fetal stages, which underpinned fetal 
myofiber hyperplasia and postnatal muscle hypertrophy. Moreover, we identified some differentially expressed 
genes between the two breeds that could be potential myostatin targets for further investigation. 

Conclusions: Proliferation of myofibers proceeded in a pulsed-wave mode at particular fetal stages in the sheep. 
The myostatin mutation changed the gene expression pattern in skeletal muscle at a transcriptome-wide level, 
resulting in variation in myofiber phenotype between Texel and Ujumqin sheep during the second half of 
gestation. Our findings provide a novel and dynamic description of the effect of myostatin on skeletal muscle 
development, which contributes to understanding the biology of muscle development in large mammals. 




Genomics 



Background 

Texel sheep, a typical "double muscle" breed due to a 
GDF8 mutation [1-3], are now commercially produced 
throughout the world, with no adverse effects detected 
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by objective assessments of meat quality [4]. However, 
evidence for an association between g+6723G > A and 
decreased intramuscular fat and reduced eating quality 
has been observed [5], Compared with Texel sheep, 
indigenous Chinese Ujumqin sheep, with no GDF8 
mutation [6], are less muscular and have a higher fat 
content, but they are superior in terms of perceived 
meat quality. Therefore, these two sheep breeds provide 
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a good natural model for studying muscle and fat devel- 
opment, as well as for identifying myostatin genes. 

Prenatal skeletal muscle development is an important 
determinant of both muscularity and meat quality [7]. 
In large precocial species such as sheep [8,9] and cattle 
[10], the maximum myofiber complement of a muscle is 
achieved prior to birth. More than three waves of myo- 
genic cells appear in sheep, and most myofibers form 
during the second half of gestation [11,12]. However, 
whether the myofibers increase with a pulsed-wave 
mode at certain developmental stages or whether they 
augment evenly across developmental stages in fetal 
sheep remains unclear. 

Myostatin, a member of the transforming growth fac- 
tor-|3 (TGF-P) family, is predominantly expressed and 
secreted by skeletal muscle and functions as a negative 
regulator of muscle growth. Mutations in the myostatin 
gene lead to a hypertrophic phenotype in mice, sheep, 
cattle, dog, and human [1,3,5,13-18]. The effect of myos- 
tatin on gene expression in prenatal muscles at the gen- 
ome-wide level was recently explored in fetal cattle 
[19-22], but no studies have been conducted dynami- 
cally at multiple fetal stages comparing two pure breeds 
with extreme phenotypes. A recent mice study demon- 
strated that myoblasts from embryonic and fetal stages 
not only had different fusion abilities, proliferation, dif- 
ferentiation and responses to TGF-P , phorbol ester 12- 
O-tetradecanoylphorbol-13-acetate, and bone morphoge- 
netic protein-4 in vitro, but they also differed in gene 
expression profiles [23], indicating that complicated and 
obvious changes in physiology and biochemistry occur 
during the prenatal stage in vivo. Therefore, investigat- 
ing the subtle changes of the effect of a myostatin muta- 
tion on skeletal muscle development at multiple fetal 
stages using our experimental model is necessary. 

Here, we examined gene expression and myofiber 
development in Texel and Ujumqin sheep at the tran- 
scriptome and histological levels. Our findings contri- 
bute to understanding the dynamic effects of myostatin 
on the biology of prenatal skeletal muscle development 
in large mammals. The results also provide clues into 
human myopathy and obesity during prenatal stages. 
We also identified putative candidate genes that are 
valuable for meat-quality traits in farm animals. 

Results 

Histological characteristics of ovine fetal skeletal muscle 
development 

We examined the number and diameter of myofibers at 
each developmental stage in Texel and Ujumqin sheep 
via histological analysis (hematoxylin and eosin [H and 
E] staining) (Figures 1 and 2). Significant difference was 
observed in myofiber diameter between Texel and 
Ujumqin sheep at each development stage, except at 70 
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Figure 1 Changes in muscle fiber diameter throughout the five 
developmental stages in Texel and Ujumqin sheep Longissimus 
dorsi muscle fiber diameter was examined in cross sections using 
hematoxylin and eosin staining. At least three animals were used to 
analyze each development stage in each breed. Photomicrographs 
of five to eight randomly selected areas were used to calculate the 
diameter of myofibers in each cross section of muscle examined. 



d. However, the myofiber diameter patterns of the two 
breeds were similar throughout the five developmental 
stages. 

Unexpectedly, we found that the myofiber hyperplasia 
profiles differed between Texel and Ujumqin fetuses 
(Figure 2). The proliferation of myofibers in Texel sheep 
was fastest at 70 and 85 d, particularly at 85 d, suggest- 
ing that this period is vital for muscle development in 
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Figure 2 Changes in muscle fiber number throughout the five 
developmental stages in Texel and Ujumqin sheep Longissimus 
dorsi muscle fiber number was examined in cross sections using 
hematoxylin and eosin staining. At least three animals were used to 
analyze each development stage in each breed. Photomicrographs 
of five to eight randomly selected areas were used to calculate the 
total number of myofibers in each cross section of muscle 
examined. 
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fetal Texel sheep. Myofiber hyperplasia reached a peak 
at 100 d in Ujumqin sheep, when the number of myofi- 
bers exceeded that of Texel sheep. Notably, with devel- 
opment, the number of myofibers in Texel fetuses 
recovered to that of Ujumqin fetuses at 120 d, although 
myofiber hypertrophy in Texel sheep was greater than 
that in Ujumqin sheep at this stage. Therefore, a new 
wave of myogenesis most probably occurs at 120 d in 
Texel sheep, given that the number of myofibers 
increased in Ujumqin sheep during this period. 

Additionally, many small-diameter fibers occurred in 
the muscle at 135 d in both breeds. The small fibers 
could be due either to intrafascicular terminations [24] 
or have been generated by the fusion of myoblasts dur- 
ing the period investigated. Therefore, to explore the 
myogenic potential in the two breeds of sheep, we iden- 
tified Pax7-positive cells in skeletal muscle at 135 d 
using immnofiuorescence staining (IFS) and found that 
the ratio of the number of Pax7-positive cells (myoblasts 
or satellite cells) to the fibers on each muscle slide was 
significantly higher in the Texel fetuses than in the 
Ujumqin fetuses (Figure 3). In combination with Figure 
2, this observation suggests that a period of 100-120 d is 
another new proliferative stage for myofiber hyperplasia 
in Texel sheep. Higher numbers of Pax7-positive cells 
were correlated with a higher muscle mass in Texel 
sheep during both the prenatal and postnatal stages. 
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Figure 3 Identification of Pax7-positive cells in longissimus 
dorsi muscle by immunofluorescence staining Longissimus 
tissue sections were stained with fluorescent antibodies Pax7 
(green) and laminin (red) at 135 d of gestation, respectively. A-D are 
the same frozen section stained with the two different fluorescent 
antibodies (400x). The corresponding antigens in muscle are 
indicated by an arrow (Pax7) and an arrowhead (laminin). 



Microarray experiment 

A mean transcript sequence identity of 97% was found 
between ovine and bovine orthologs, which probably 
highlights the utility of the microarray technique with 
ovine samples [25]. However, some no-reporting probe 
sets were found, partly because of the poor performance 
of the probe sets with samples from a related species. 
The best way to resolve this type of problem is to use a 
species-specific microarray. In the present study, we 
applied the first specialized and standardized transcrip- 
tome-wide sheep oligo DNA microarray (Agilent Sheep 
Gene Expression Microarray; Agilent Technologies, 
Santa Clara, CA, USA). Due to improvements in species 
specificity, the average detectable rate of all probe sets 
reached 87.83% (Additional file 1), which enabled us to 
capture the subtle changes in gene expression. 

A pool of differentially expressed (DE) probes from 
the two groups was used for a systematic hierarchical 
clustering to gain insight into the transcriptome-wide 
similarities among all 31 individuals investigated (Figure 
4). We found that individuals at the same developmental 




Figure 4 Hierarchical clustering between the differentially 
expressed genes and individuals from the two sheep breeds 

Clustering was performed using GeneSpring 10.0. A one-way 
analysis of variance was applied to two within-breed contrasts 
across developmental time. A pool of the differentially expressed 
probes from the two groups was used for system hierarchical 
clustering to clarify the transcriptome-wide similarities among all 31 
individuals investigated. 
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stage were clustered together regardless of their genetic 
background, indicating that differences derived during 
development were larger than those from the genetic 
backgrounds between the breeds. 

Gene expression in skeletal muscle at various 
developmental stages between breeds 

According to the primary functions of different genes 
including muscle, lipid, and the immune system, we 
visualized the data on MA and bar plots (Figures 5 and 
6). Combined with Figures 1 and 2, Figure 6 also sug- 
gests the phenotypic profile in muscle and the profile of 
a given category of DE genes at different developmental 
stages. The biological functions and canonical pathways 
were explored at various developmental stages in Texel 
and Ujumqin sheep (Tables 1, 2, 3, 4 and 5). The partial 
DE genes in the two breeds of sheep at various develop- 
mental stages are listed in Tables 6 and 7 (the full set of 
differential genes at each developmental stage between 
the breeds is listed in Additional files 2, 3, 4, 5 and 6). 
170 vs. U70 

The rapid increase in total muscle fiber number in the 
ovine fetus begins at 70 d [11,12]. We found that the 
number of DE genes was highest at this stage compared 



with other stages: 207 probe sets, including 155 known 
genes and 50 expressed sequence tags (ESTs) were dif- 
ferentially expressed in the T70 and U70 samples. Com- 
pared with U70, 82 genes were upregulated and 73 
genes were downregulated in T70. A gene ontology ana- 
lysis by DAVID revealed that these upregulated genes 
encoded cell activation, negative regulation of apoptosis, 
positive regulation of peptidyl-tyrosine phosphorylation, 
and vesicle-mediated transport. Genes that were down- 
regulated encoded hexose catabolism, monosaccharide 
catabolism, postembryonic development, alcohol catabo- 
lism, and lipid biosynthesis. 

To ascertain whether some molecular events or cas- 
cades were associated with the process of skeletal mus- 
cle development, we identified four canonical pathways 
by DAVID and the IPA platforms, which are involved in 
myogenesis and muscle metabolism, including the glyco- 
lysis/gluconeogenesis pathway, Af-glycan biosynthesis, 
calcium signaling, and chemokine (C-X-C motif) recep- 
tor 4 (CXCR4) signaling). 

We also revealed through IPA software that genes 
associated with the skeletal and muscular system devel- 
opment and function, as well as lipid metabolism, were 
top biological functions (Table 1). 
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Figure 5 MA plot for differential genes in Texel and Ujumqin sheep. Genes expressed more highly in Ujumqin sheep are on the bottom, 
and genes expressed more highly in Texel sheep are on top. Genes involved in muscles, lipids, and the immune system are denoted by blue 
triangles, green circles, and red rectangles, respectively. For example, HLA-DQB1 and AR were highly expressed in Texel, whereas MYL6B was 
highly expressed in Ujumqin sheep. Other genes expressed differentially are denoted by gray rectangles. 
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Figure 6 Abundance of primary functional clustering of differentially expressed genes Counts of genes involved in muscles, lipids, and 
the immune system are represented by the ratio of functional categories relative to the total number of genes that exists in the microarray chip 
for each of those categories, respectively. The ratio illustrates the abundance of different clusters of differentially expressed genes between Texel 
and Ujumqin sheep at five developmental stages (70, 85, 100, 120, and 135 d of gestation). 



We further analyzed some transcripts and found that 
eukaryotic translation elongation factor 1 alpha 2 
(EEF1A2) was highly expressed in T70. In contrast, the 
expression levels of signal transducer and activator of 
transcription 3 {STATS), aldolase A (ALDOA), and the 
leptin receptor {LEPR) were upregulated in U70. 
785 vs. U85 

In total, 115 probe sets were differentially expressed in 
T85 and U85, including 97 annotated genes and 18 
ESTs. In contrast to U85, 50 genes were upregulated 
and 47 were downregulated in T85. Of these genes, 
immune response, cellular carbohydrate catabolism, and 
fatty acid metabolism were significantly overrepresented 
in the upregulated genes, whereas cell migration regula- 
tion was significantly overrepresented in the downregu- 
lated genes. 

For all DE genes between U85 and T85, we identified 
four canonical pathways that were overrepresented 
between the two breeds: glycolysis/gluconeogenesis, 
fructose and mannose metabolism, the pentose phos- 
phate pathway, and vascular endothelial growth factor 
(VEGF) signaling. Moreover, hematological system 
development and function and lipid metabolism were 
still significantly represented as top biofunctions. In par- 
ticular, we found a cluster of genes enriched in the bio- 
function of skeletal and muscular disorders, suggesting a 
difference in skeletal muscle development between 
Texel and Ujumqin fetal sheep (Table 2). 

Among the DE genes, eukaryotic translation elonga- 
tion factor 1 alpha 2 {EEF1A2), fructose-1,6- 



bisphosphatase 2 {FBP2), actinin alpha 3, proteoglycan 3 
{PRG3), myosin heavy chain 1 {MYH1) were more 
highly expressed in Texel than in Ujumqin sheep. How- 
ever, the stathmin-like 2 {STMN2) and phosphoinosi- 
tide-3-kinase, regulatory subunit 1 alpha {PIK3R1) were 
more highly expressed in U85. 
7700 vs. U100 

At this fetal developmental stage, the difference in mus- 
cle development, particularly the rate of myofiber prolif- 
eration, was most significant between the breeds. Ninety 
probe sets, including 70 annotated genes and 20 ESTs, 
were differentially expressed in T100 and U100. Among 
these genes, 33 were upregulated and 37 were downre- 
gulated in T100. The upregulated genes were involved 
in DNA binding, RNA binding, and ribonucleoproteins, 
whereas those that were downregulated were mostly 
related to steroid biosynthesis, lipid metabolism, and 
skeletal muscle fiber development. 

For all DE genes at 100 d, hematological system devel- 
opment and function remained the most overrepre- 
sented physiological system, and lipid metabolism was 
also more active at this stage. Inositol metabolism, ster- 
oid biosynthesis, eukaryotic initiation factor 2 (EIF2) sig- 
naling, and VEGF signaling were overrepresented 
between T100 and U100 (Table 3). 

Within the DE genes, thyroid stimulating hormone 
receptor {TSHR) and paraoxonase 1 {PON1) were more 
highly expressed in Ujumqin sheep than in Texel sheep, 
whereas OVAR-DQB1 (sheep MHC class II), pyruvate 
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Table 1 Summary of the differentially expressed genes between T70 and U70 using IPA 

Top Bio Functions 
Diseases and Disorders 



Name 


p-value 1 


Molecules 2 
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51 
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19 


Developmental Disorder 
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R 91 F-06-7 93F-03 


8 


nflammatory Response 


1 .05E-05-7.93E-03 


25 


Molecular and Cellular Functions 






Name 


p-value 
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Cpllular Mnupmpnt 
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Lipid Metabolism 


9.16E-06-7.93E-03 


25 
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28 


Small Molecule Biochemistry 


9.16E-06-7.93E-03 


30 


Physiological System Development and Function 






Name 


p-value 


Molecules 


Tissue Development 


1 .77E-06-7.93E-03 


24 


Hematological System Development and Function 


7.13E-06-7.93E-03 


22 


Hematopoiesis 


8.91 E-06-7.93E-03 


9 


Immune Cell Trafficking 


8.91 E-06-7.93E-03 


19 


Skeletal and Muscular System Development and Function 


2.68E-05-7.93E-03 


18 


Top Canonical Pathways 






Name 


p-value 


Ratio 3 


Acute Phase Response Signaling 


2.2E-03 


6/1 78 (0.034) 


Glycolysis/Gluconeogenesis 


7.65 E-03 


4/142 (0.028) 


A/-Glycan Biosynthesis 


9.44E-03 


3/93 (0.032) 


Calcium Signaling 


1.17E-02 


5/204 (0.025) 


CXCR4 Signaling 


3.31-02 


4/1 69 (0.024) 



1. The p-value was calculated using a right-tailed Fisher's exact test. 2. The number of genes enriched in corresponding functions. 3. The ratio was calculated as 
the number of genes from our data set that overlap with the canonical pathway in question divided by the total number of genes represented in that canonical 
pathway. 

4. Italicized terms such as "Cancer" and "Cardiovascular Disease" represent a biased annotation or are not related to the aim of the study. 
IPA version: 8.7 (release date: 09-17-2010) 
Content version: 3202 (release date: 07-07-2010) 
Reference set: Ingenuity knowledge base (genes only} 
Relationships to include: direct and indirect 
Includes endogenous chemicals 

Filter summary: Consider all molecules and/or relationships 

dehydrogenase kinase, isozyme 4 (PDK4), and thiore- 
doxin reductase (TXNRD1) were upregulated in T100. 
T120 vs. U120 

The discrepancy in fetal muscle development was more 
significant between Texel and Ujumqin sheep at 120 d. 
In total, 151 probe sets, including 120 annotated genes 
and 31 ESTs, were differentially expressed in T120 and 
U120. Of these genes, 63 were upregulated in T120 and 
57 were upregulated in U120. The upregulated genes 
were specific to the immune response, leukocyte activa- 
tion, and the I-kappaB kinase/NF-kappaB cascade, 
whereas the downregulated genes were significantly 
associated with the response to ethanol, the enzyme- 



linked receptor protein signaling pathway, and histone 
H2A acetylation. 

Hematological system development and function, 
immune cell trafficking, and nervous system develop- 
ment and function were most overrepresented in phy- 
siological systems development and function. 
Furthermore, lipid metabolism was also the most over- 
represented molecular and cellular function. Several 
canonical pathways were associated with immune and 
lipid metabolism disorders, such as nuclear factor (NF)- 
kB signaling and atherosclerosis signaling (Table 4). 

Among DE genes between T120 and U120, OVAR- 
DQB1 (sheep MHC class II), guanylate-binding protein 
2; interferon-inducible, guanylate-binding protein 4; and 
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Table 2 Summary of the differentially expressed genes between T85 and U85 using IPA 

Top Bio Functions 
Diseases and Disorders 
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Endocrine System Development and Function 


3.71E-04-1.98E-02 


5 


Visual System Development and Function 
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Hematopoiesis 
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8 


Immune Cell Trafficking 
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16 


Top Canonical Pathways 






Name 


p-value 


Ratio 3 


Glycolysis/Gluconeogenesis 


1 34E-04 


5/142 (0.035) 


Fructose and Mannose Metabolism 


2.19E-04 


4/145 (0.028) 


Pentose Phosphate Pathway 


9.13E-04 


3/89 (0.034) 


Dendritic Cell Maturation 


1.31E-03 


5/188 (0.027) 


VEGF Signaling 


9.55E-03 


3/99 (0.03) 



Symbols and analytical settings are the same as in Table 1. 



CCAAT/enhancer-binding protein (C/EBP) beta 
{CEBPB) were upregulated in T120. Both STMN2 and 
transformation/transcription domain-associated protein 
(TRRAP) were upregulated in U120. 
T135 vs. U135 

Seventy-two probe sets, including 58 annotated genes 
and 14 ESTs, were differentially expressed in T135 and 
U135. Of these genes, 27 were upregulated in T135 and 
another 31 were more highly expressed in U135. The 
immune response process was overrepresented in the 
upregulated genes, and mitosis was represented signifi- 
cantly in the downregulated genes. 

For the pool of DE genes, embryonic development and 
tissue development, hematological system development 
and function, cell-mediated immune response, and ner- 
vous system development and function were listed at 
the top of physiological system development and func- 
tion. Cellular functions such as the cell cycle, cell death, 
and cell-to-cell signaling and interaction were more 
active. Notably, a pathway involved in glycerolipid meta- 
bolism was identified between T135 and U135 (Table 



5), two members of which (ALDH1A1 and APOC2) 
were upregulated in U135. 

In terms of physiological system development and 
functions, differences in the immune and hematological 
systems were most significant between Texlel and 
Ujumqin sheep during muscle development. Then, dif- 
ferences in nervous system development and function 
became more prominent at later gestational stages. The 
most obvious differences in molecular and cellular func- 
tions were lipid metabolism and cell communication 
between the two breeds. Several overrepresented canoni- 
cal pathways related to muscle and adipose development 
(Table 5), which were particularly noteworthy, under- 
pinned the differences in skeletal muscle and adipose 
tissue between the two breeds at the fetal stage. 

Quantitative PCR validation of microarray gene 
expression 

To confirm the differential expression of the genes in 
the microarray analysis, we selected 15 genes to validate 
by quantitative PCR (qPCR). Among these, eight genes 
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Table 3 Summary of differentially expressed genes between T100 and U100 using IPA 
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7 08F-04-3 S7F-0? 

J.J/ l_ \JjL 


] ] 




5.03E-04-1.92E-02 


M 


fVlnlpn Mar anH ("pllnlar Funrtinnc 






Name 


p-value 


Molecules 


Lipid Metabolism 


2.91 E-1 1-3.57E-02 


16 


Sma Molecule Biochemistry 


2.91 E-1 1-3.57E-02 


21 


Vitamin and Mineral Metabolism 


2.91 E-1 1-2.60E-02 


16 


Mnlpriil^r Tr.^nc.nnrt 

i V 1 U 1 Q 1 1 1 a 1 1 J LJU 1 L 


8 41F-06-3 S7F-02 

Cj.*-r 1 l_ UU J.J/ l_ 


1 9 


C pIIi liar (\/lru/pmpnt 
v^CI 1 Ulcil iviUvtriiltrML 


Z.UOL W J ) / |_ KJA. 


5 


PKucinlnn ira 1 Cuctom nouolnnmont anH Fi mrt'mn 

rllyjIUIUUILCII jySLCNI L/C VCIULJI 1 ICI 1 1 aiiu IUI ILIIUI 1 






Name 


p-value 


Molecules 


Hpmstnlnn irs 1 ^v^tprn Flpvplnnmpnt ^nH Fi inrtinn 

1 ICTIIIQLUILy^JI^Gll J \0 LCT 1 1 1 LJtTVCIVJlJIIItrlll a 1 1 U 1 U 1 lv~LI\JI 1 


Z..L^CJl_ J.J/ l_ \JjL 


1 1 


Immune Cell Trafficking 


2.08E-04-3.57E-02 


9 


Lymphoid Tissue Structure and Development 


2.08E-04-3.57E-02 


4 


Tissue Development 


479E-04-3.57E-02 


9 


Hematopoiesis 


6.96E-04-3.57E-02 


6 


Top Canonical Pathways 






Name 


p-value 


Ratio 3 


T-Helper Cell Differentiation 


2E-03 


3/72 (0.042) 


Inositol Metabolism 


373E-03 


3/97 (0.031) 


Biosynthesis of Steroids 


5.98E-03 


2/128 (0.016) 


EIF2 Signaling 


3.89E-02 


2/103(0.019) 


VEGF Signaling 


3.98E-02 


2/99 (0.02) 



Symbols and analytical settings are the same as in Table 1. 



(ACTB, TNC, MY09A, MYL6B, MYH1, CASQ1, TPM2, 
and TMOD4) encode myofibrillar proteins and two 
(EEF1A2 and LEF-1) are transcription factors that regu- 
late the genes associated with myogenesis. PIK3R1 is 
involved in fat metabolism regulation, whereas TRRAP, 
TXNRD1, and INSIG1 participate in the regulation of 
gene expression and cell proliferation. Additionally, we 
also detected MSTN expression during skeletal muscle 
development. As in previous studies [26-29], RpLPO was 
chosen as an ideal reference gene to normalize the data. 
The ratios of 12 DE genes at particular developmental 
stages in the two sheep breeds in the microarray analysis 
agreed with the qPCR results (Table 8). 

Discussion 

Characteristics of prenatal muscle fiber development in 
sheep 

Wilson et al. and Maier et al. showed that at least three 
waves of myogenesis occurred until d 76 of gestation in 
sheep [9,11]. We demonstrated that the myofiber num- 
bers increase sharply with a pulsed-wave mode at a 



particular developmental stage in sheep. This is the first 
time that the mode of fetal myofiber proliferation was 
revealed in a large mammal during the second half of 
gestation. In addition to the longissimus dorsi muscle, 
similar patterns of myofiber proliferation were found in 
the semitendenosus, gluteus medius, gastrocnemius, and 
triceps brachii muscles in Texel and Ujumqin fetal 
sheep (data not shown). These findings indicate that the 
specific developmental stage during which myogenic 
surges of myofiber hyperplasia appeared is pivotal for 
fetal muscle development as well as postnatal muscular- 
ity potential in sheep. 

Unlike the studies of Wilson et al. and Maier et al., 
who identified different generations of myogenic cells 
using electron microscopy [9,11], we counted the muscle 
fibers under a light microscope and found many small- 
diameter fibers until 135 d in both breeds. We consid- 
ered that at least part of the small fibers must have 
been generated by the fusion of myoblasts (or satellite 
cells) during the study period. Pax7 plays an important 
role conferring myogenic potential to these progenitor 
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Table 4 Summary of differentially expressed genes between T120 and U120 using IPA 

Top Bio Functions 
Diseases and Disorders 



Name 


p-value 1 




Molecules 2 


Hypersensitivity Response 


9.54E-06-1 


19E-02 


7 


Inflsmmsfnrv Rp^nnn^p 

II II 1 □ 1 1 II 1 laLUI y IiCjLJWI IjCT 


9.54E-06-1 


1 9E-02 


18 


Cancer 4 


3.78E-05-1 


.1 9E-02 


41 


1— Ipmstnlnn irs 1 Pli^p^^p 

1 Id 1 1 a LU 1 KJ\^ 1 t-Q 1 L/ljCra_iC: 


S 66F-0S-1 

J.UUL WJ 1 


1 9E-02 


7 


Immunologica Disease 


5.66E-05-1 


.1 9E-02 


7 


Mnlpn ilar anrl ("pllnlar Funrtinnc 

iviuicluiui ci i iu uciiuiui rui 1 l 1 KJ 1 I J 








Name 


p-value 




Molecules 


Cell - To - Cell Signaling and Interaction 


9.54E-06-1 


.1 9E-02 


17 


Lipid Metabolism 


1 .1 6E-05-1 


.1 9E-02 


16 


Small Molecule Biochemistry 


1 .1 6E-05-1 


.1 9E-02 


20 


rai"hnh\/n rpfp K/lpfaKn Km 
v_di uui i y u i ate i vitriaLJiJi i ii i i 


3 54E-05-1 


1 9E-02 


5 


Dn in MptnhnlKm 


3.54E-05-1 


1 9E-02 


1 0 


Phu^inlnniral ^v^tpm Rpvplnnmpnt and Funrtinn 

r 1 IV JlUIUUIVOI JyjlUI 1 1 L/C VCIUUIIICIIl ul IU 1 Ul ILUUI 1 








Name 


n-v^ 1 1 ip 
yj value 




Mnlpn lips 

IVIUIC^UICj 


Hematological System Development and Function 


9.54E-06-1 


19E-02 


23 


Immune Cell Trafficking 


9.54E-06-1 


19E-02 


13 


Nervous System Development and Function 


1.69E-05-1 


.19E-02 


8 


Tissue Development 


1.69E-05-1 


.19E-02 


13 


Reproductive System Development and Function 4 


7.17E-05-1 


.19E-02 


10 


Top Canonical Pathways 








Name 


p-value 




Ratio 3 


Altered T Cell and B Ceil Signaling in Rheumatoid Arthritis 4 


1 .54E-04 




5/91 (0.055) 


NF-kB Signaling 


1 .87E-03 




5/155 (0.032) 


Autoimmune Thyroid Disease Signaling 


2.77E-03 




3/61 (0.049) 


Atherosclerosis Signaling 4 


2.92E-03 




4/113 (0.035) 


Taurine and Hypotaurine Metabolism 


3.53E-03 




2/47 (0.043) 



Symbols and analytical settings are the same as in Table 1. 



cells [30], and Pax 7 is an ideal cell marker of myogenic 
progenitor cells such as myoblasts and satellite cells 
[31]. Our findings indicate a higher potential for postna- 
tal muscle growth that entails more Pax7-positive cells 
in the fetal stage. 

Gene expression in the two breeds at various 
developmental stages 

Of the five prenatal stages, the numbers of DE genes in 
the two sheep breeds were highest at 70 and 120 d, indi- 
cating that rapid myofiber proliferation included numer- 
ous genes during the myogenic process. However, two 
surges of myofiber proliferation appeared at 85 and 120 
d in Texel fetuses, probably because the expression of 
more genes associated with myoblast proliferation and 
differentiation at 70 d prepared for the myofiber hyper- 
plasia that followed at 85 d, and more genes related to 
the immune system were differentially expressed 
between T70 and U70 (Figure 6), which led to larger 
total counts of DE genes. The difference between the 
two sheep breeds in the timing of myogenesis in terms 



of transcriptomic levels was similar to that found 
between pig breeds [32]. Cagnazzo et al. showed that 
myogenesis-related gene expression is greater in early 
Duroc (a breed with more intramuscular fat) embryos 
than in early Pietrain (a highly muscled breed) embryos 
at 14-49 d of gestation, whereas the opposite was found 
in late embryos (63-91 d of gestation) [32]. The present 
findings suggest that highly muscled breeds have a 
longer myogenic process during prenatal stages and that 
myogenesis is more intense in late-stage fetuses, which 
was validated by Pax7-cell staining in longissimus mus- 
cle cross sections (Figure 3). However, whether this 
myogenic process is the same in sheep as in pigs at ear- 
lier embryonic stages remains to be investigated. 

We found that the immune and hematological sys- 
tems' development and function were most overrepre- 
sented of all physiological systems and functions, 
accompanied by muscle system development, which was 
consistent with a microarray developmental analysis in 
cattle fetuses [21]. Recent discoveries have revealed 
complex interactions between skeletal muscle and the 
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Table 5 Summary of differentially expressed genes between T135 and U135 using IPA 

Top Bio Functions 
Diseases and Disorders 



Name 


p-value 1 


Molecules 2 


Cancer^ 


9.23E-06-2.87E-02 


7 


Rpnrnrlt irth/p 'sv^tpm Dl^pn^p^ 


9 73F-06-6 44F-03 

J.Z.JL \j\J vJ.T^L UJ 


7 


Genetic Disorder 


3.22E-03-2.87E-02 


5 


mm innlnnir^l PIkp^p 

II 1 II \ 1 U 1 IUIUU ILal \-J 1 jCTu jCT 


3 22F-03-3 22F-03 


1 


nfection Mechanism 


3.22E-03-4.1 1 E-02 


4 


Mnlpn Mar and ("pllular Funrtinnc 






Name 


p-value 


Molecules 


Cell Cycle 


5.50E-04-2.87E-02 


3 


r pii np^rh 

\— Cri I \-J\ZO. LI 1 


1 04F-03-? 43F-07 


8 


Cell-To-Cell Signaling and Interaction 


3.22E-03-4.42E-02 


5 


Cpllulsr A^^pmhlu ^nH Orn^ini73tirin 

\— Cri 1 LA 1 Q 1 mjjCI 1 l\*Jl y O 1 1 LJ W ivjai l lz_u i_ 1 W 1 1 


3 77F-03-442F-0? 

-J .Z.Z L_ WJ ^ r^Z. I_ W Z 


7 


C p i ibr C nmnrnmicp 
v^CI 1 UICII v_Ul 1 IkJI \J\ 1 II JC 


3 99F-n^-44.?F-D9 

J./ZL UJ TVTiL. UZ 


3 


Ph\/cifilnn \C7\ 1 ^mtpm Dp\/pIo nm ont anH Fi inrtirtn 

rllyjIUIUUILCII jySLCNI L/C VCIUIJI 1 Id 1 I allU lUIILLIUII 






Name 


n-valiiP 

L/ VulUC 


Mnlprulpc 


Fmhrvnnir r~)p\/plnnmpnf 

L_IIILJiyLJIIIl~ L>tTVtriV^JLJIIIC:IIL 


5 50F-04-? SSF-02 

J.JWL Z. J JL WZ 


5 


Tissue Development 


5.50E-04-4.1 1 E-02 


7 


Hematological System Development and Function 


247E-03-4.42E-02 


7 


Cell-mediated Immune Response 


322E-03-4.42E-02 


2 


Nervous System Development and Function 


322E-03-6.44E-03 


2 


Top Canonical Pathways 






Name 


p-value 


Ratio 3 


LPS/IL-1 -mediated Inhibition of RXR Function 


3.61 E-03 


4/216 (0.019) 


Histidine Metabolism 


1 26E-02 


2/120 (0.017) 


PXR/RXR Activation 


22E-02 


2/91 (0.022) 


Arginine and Proline Metabolism 


2.68E-02 


2/183 (0.011) 


Glycerolipid Metabolism 


4.51 E-02 


2/156 (0.013) 



Symbols and analytical settings are the same as in Table 1. 



immune system that regulate muscle regeneration and 
myogenesis, and many immune molecules, such as 
tumor necrosis factor-alpha (TNF-a), NF-kB, interleukin 
(IL)-4, IL-6, IL-10, and leukemia inhibitory factor are 
involved in muscle cell proliferation and differentiation 
[33-38]. In most cases, whether perturbation in the 
same signaling pathways related to myogenesis occurs 
will determine normal myogenic development or muscle 
disorders. Therefore, the pathways associated with 
immune and muscle disorders in this work are valuable 
for human myopathy at the prenatal stage. The hemato- 
logical system is also implicated in myogenesis [39-43], 
which probably reflects the systemic requirement for 
muscle function during fetal development. 

Nervous system development was overrepresented in 
the last two development stages in the two sheep breeds. 
Neurons are indispensable for maintaining normal mus- 
cle physiological function and also parallel the skeletal 
muscle in auxology. Comprehensive and close interac- 
tions occur between muscles and the nervous system 
during development [44-46]. Muscle development is 



regulated by the central nervous system in Drosophila 
and pigs [47,48], and muscle fiber type is dependent on 
the pattern of innervation of a muscle established due 
to differential projection patterns between fast and slow 
motoneurons [49,50]. Animals with different muscle 
phenotypes undergo diverse innervation patterns during 
fetal development. Double-muscled cattle have an addi- 
tional 13-26% increased branching in terminal axons 
compared to that in normal cattle caused by a real 
increase in the number of myofibers [51]. 

We identified several valuable canonical pathways 
directly associated with muscle development and func- 
tion, myogenesis, myoblast proliferation, and the cell 
cycle in muscle at various developmental stages, such as 
calcium [52], CXCR4 [53-55], and VEGF signaling [56]. 
Furthermore, other pathways involved in adipose and 
muscle metabolism, such as inositol metabolism, steroid 
biosynthesis, EIF2 signaling, glycerolipid metabolism, the 
glycolysis/gluconeogenesis pathway, and AT-glycan bio- 
synthesis, are significantly affected during development. 
Steelman et al. suggested that Wnt signaling is a 
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Table 6 Examples of downregulated genes between Texel and Ujumqin sheep at the same developmental stage 



Probe Name 


Genbank Accession 
Number 


Gene 


rolu 
Change 


Developmental 
stage 


A. 


70. 


.P069356 


EE752033 


Homo sapiens selenoprotein I (SELI) 


15.35 


70 d 


A 

rt_ 


7fl 
_/ U_ 


_ru i ozy I 


EE793140 


Homo sapiens Rho-related BTB domain containing 3 (RHOBTB3) 


5.16 


70 d 


/\ 


70 


POOASAI 
_r uuuju i 


EE781286 


Homo sapiens nucleolar protein 7 (NOL7) 


4.54 


70 d 


A 

A_ 


_/u_ 


_rUU'HDZ I 


EE863207 


Homo sapiens signal transducer and activator of transcription 3 (STAT3) 


2.38 


70 d 


A 

A_ 


_/u_ 


Dnm 7H£ 
_rUU I /UD 


FE033261 


Homo sapiens aldolase A, fructose-bisphosphate, transcript variant 1 (ALDOA) 


2.28 


70 d 


A 


7fl 

_/ U_ 


jUjOo 1 Z 


NM_001 009763 


Homo sapiens leptin receptor (LEPR) 


2.24 


70 d 




70 


_r uujjzu 


CO202686 


Homo sapiens stathmin-like 2 (STMN2) 


8.23 


85 d 


A 


7n 


_ru i ozy i 


EE793140 


Homo sapiens Rho-related BTB domain containing 3 (RHOBTB3) 


5.99 


85 d 


^ 


70 


juoo i yo 


EE779349 


Homo sapiens trefoil factor 3 (TFF3) 


5.1 9 


85 d 


A 

AA_ 


7n 


_rUjyjj 1 


NM_00 1009266 


Homo sapiens proopiomelanocortin (POMC) 


5.00 


85 d 


A 

A_ 


7n 


_rUj /uyo 


EE773353 


it ■ u 1" * . ■ i i l ■ I4- 1 '.1/11 \ /nii/ini \ 

Homo sapiens phosphoinositide-3-kinase, regulatory subunit 1 (alpha) (PIK3R1), 
transcript variant 2 


2.77 


85 d 


A 


_/ u_ 


pnni 731 


EE780294 


Homo sapiens trafficking protein, kinesin-binding 1 (TRAK1) 


8.04 


100 d 


A 

r\_ 


7n 
_/u_ 


DH1 701 A 
_rU I / y I D 


EE828626 


Homo sapiens S100 calcium-binding protein Al 1 (SI 00A1 1 ) 


8.03 


100 d 


A 


70 


P032741 


CF1 16713 


Bos taurus derensin, beta (DEFB) 


7.09 


100 d 




70 


P00141 1 


EE873259 


805 taurus secretoglobin, family 1D, member 2 (SCGB1D2) 


5.59 


100 d 




70 


P039326 


NM_00 1009410 


Homo sapiens thyroid-stimulating hormone receptor (TSHR) 


4.50 


100 d 




70 


P045271 


r~ 1 — 7nm f\r\ 

EE788300 


Homo sapiens paraoxonase 1 (PON1) 


4.01 


100 d 


A 

A_ 


7n 
J u_ 


dh^q 1 3 1 

Y UOo I J I 


DQ1 52973 


Homo sapiens lymphoid enhancer-binding factor 1 (LEF-1) 


2.23 


100 d 


A. 


70. 


.P035741 


FE033512 


Bos taurus X (inactive)-speciric transcript (XIST), noncoding RNA 


65.59 


120 d 


A. 


70. 


.P005526 


CO202686 


Homo sapiens stathmin-like 2 (STMN2) 


21.27 


120 d 


A. 


70. 


_P0 18856 


EE791277 


Homo sapiens ATP-binding cassette, subfamily B (MDR/fAP), member 4 
(ABCB4) 


15.37 


120 d 


A. 


70. 


.P02481 1 


EE777383 


Homo sapiens epithelial cell adhesion molecule (EPCAM) 


7.1 6 


1 20 d 


A. 


70. 


_P0 19481 


FF89881 3 

LLOZOO I J 


\-ir\rnr\ (/in/on c trancfnrnn^tirin/tranc^rintinn HrirTiain-3c:c(^/~iatciH nrntpin fTRRAP^ 
( i\jt i lu jup/ci fj li oi 1 51 KJi 1 1 la liu 1 1/ u dl i ji_i iuliuii i_iui l la 1 1 i a z) jULid icu uiu Lei ii \\ nnrAr j 


492 


1 20 d 


A. 


70. 


.P030147 


NM_001 009261 


Homo sapiens prion protein 2 (doublet) (PRND) 


4.84 


120 d 


A. 


70. 


.P025536 


FE022583 


Homo sapiens glutamate receptor, metabotropic 8 (GRM8) 


4.62 


120 d 


A. 


70. 


.P051317 


NM_001 009479 


Homo sapiens natriuretic peptide precursor C (NPPC) 


4.46 


120 d 


A. 


70. 


.P030402 


NM_001 093787 


Homo sapiens kallikrein 1 human (KLK1) 


4.37 


120 d 


A. 


70. 


.P051037 


NM_001 009451 


Homo sapiens solute carrier family 2 (facilitated glucose/fructose transporter), 
member 5 (SLC2A5) 


8.14 


135 d 


A. 


70. 


.P035536 


DY520412 


Homo sapiens ankyrin repeat domain 1 (cardiac muscle) (ANKRD1) 


5.37 


135 d 


A. 


70. 


_P0 16291 


EE793140 


Homo sapiens Rho-related BTB domain containing 3 (RHOBTB3) 


4.14 


135 d 



Differentially expressed probes were detected from five f-test analyses using GeneSpring 10.0 (Agilent Technologies). 

1. Original probe name in microarray analysis. 2. GenBank Accession Number of expressed sequence tags provided by the manufacturer according to the specific 
probe name. 3. The highly homologous human sequences acquired using BLASTN search in NCBI (search setting: query coverage not less than 50%, E-value less 
than 1.00E-100). 



potential downstream target of myostatin for postnatal 
skeletal muscle growth and hypertrophy in mice [57]. 
However, we focused on the transcriptome during pre- 
natal stages rather than postnatal stages and found that 
various canonical pathways are predominant in skeletal 
muscle development at different stages due to a myosta- 
tin mutation in sheep. 

DE genes involved in muscle and adipose development 

STAT3, a member of the STAT family and a cooperator 
in the Janus kinase pathway, plays a dual role in the reg- 
ulation of myoblast proliferation and differentiation, but 
is dependent on interactions with various cofactors 
[58-62]. We found that STAT3 expression was higher in 



U70 than T70, which suggests that STAT3 is most prob- 
ably associated with the differentiation rather than pro- 
liferation of myoblasts. The GO analysis of DE genes 
between T70 and U70 confirmed our presumption, and 
muscle fiber hyperplasia was more intense in T70 (Fig- 
ure 2) as well. 

ALDOA is responsible for significant activation during 
the differentiation of chicken primary myoblasts and 
plays an important role in muscle gene transcription 
[63,64]. In our study, a higher ALDOA expression level 
in U70 indicated that more myoblasts exited the cell 
cycle and entered differentiation in U70 than those in 
T70. The leptin receptor (LEPR), a protein secreted 
from adipocytes, is responsible for fat mass regulation 
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Table 7 Examples of upregulated genes between Texel and Ujumqin sheep at the same developmental stage 



Probe Name 1 


GenBank 2 Accession 
Number 


Gene^ 


roia 
Change 


Developmental 
Stage 


A. 


70. 


.P046666 


EE795873 


Homo sapiens eukaryotic translation elongation factor 1 alpha 2 (EEF1A2) 


8.92 


70 d 


A. 


70. 


.P0 11966 


FE031658 


Homo sapiens transmembrane channel-like protein 5 (TMC5) 


6.74 


70 d 


A. 


70. 


.P019721 


EE787688 


Homo sapiens complement component 4-binding protein, alpha (C4BPA) 


6.64 


70 d 


A 

r\_ 


7(1 
_/ U_ 




EE825499 


Homo sapiens CD8a molecule, complete cds (CD8A) 


6.06 


70 d 


A 


70 


_rUUj i tj 1 


CU652121 


Homo sapiens cardiomyopathy-associated 5 (CMYA5) 


4.09 


70 d 


A 


70 


_rU i +OUOO 


EE795873 


Homo sapiens eukaryotic translation elongation factor 1 alpha 2 (EEF1A2) 


16.21 


85 d 




70 




f~ 1 7~7 1 70 A 

EE775784 


Homo sapiens fructose-1,6-bisphosphatase 2 (FBP2) 


10.39 


85 d 


A 


J u_ 


PD37371 


EE776937 


Homo sapiens actinin, alpha 3 (ACTN3) 


7.83 


85 d 




70 


P054756 


EE791 183 


Homo sapiens apohpoprotein C-ll (APOC2) 


6.29 


85 d 


A. 


,70_ 


.P055056 


EE768552 


Homo sapiens proteoglycan 3 (PRG3) 


6.03 


85 d 


A_ 


,70_ 


.P050846 


LU651 71 3 


Homo sapiens myosin, heavy chain 1, skeletal muscle, adult (MYH1) 


5.60 


85 d 


A_ 


_70_ 


.P019721 


EE787688 


Homo sapiens complement component 4-binding protein, alpha (C4BPA) 


4.00 


85 d 


A_ 


.70. 


.P0 17031 


L08792 


Homo sapiens major histocompatibility complex, class II, DQ beta 1 (HLA- 
DQB1) 


131.77 


100 d 


A 
r\_ 


7n 
_/ u_ 


_r WtoU I 


EE862627 


Homo sapiens chemokine hgand 26-like protein (CCL26) 


10.16 


100 d 




70 


P024246 


tE859238 


Homo sapiens pyruvate dehydrogenase kinase, isozyme 4 (PDK4) 


5.85 


100 d 




70 


P057421 


EE784201 


Homo sapiens MAK16 homolog (S. cerevisiae) (MAK16) 


5.22 


100 d 




70 


PD7D7^A 


EE824836 


Homo sapiens thioredoxin reductase mRNA (TXNRD1) 


5.20 


100 d 


A 


70 


P01 1 371 


EE813532 


Homo sapiens isopentenyl-diphosphate delta isomerase 1 (IDI1) 


4.25 


100 d 


A 

r\_ 


J u_ 


pni 70^1 

_rU I / Uj I 


L08792 


Homo sapiens major histocompatibility complex, class II, DQ beta 1 (HLA- 
DQB1) 


1 ion 1 /" 

1389.16 


1 20 d 


A. 


.70. 


.P06601 1 


FE028749 


Homo sapiens ubiquitin-specific peptidase 9, Y- inked (USP9Y) 


67044 


120 d 


A. 


.70. 


.P066386 


EE851120 


Homo sapiens guanylate-binding protein 2, interferon-inducible (GBP2) 


76.81 


120 d 


A. 


.70. 


.P052621 


EE7781 1 1 


Homo sapiens pyrophosphatase (inorganic) 1 (PPA1) 


14.06 


120 d 




70 


P065251 


EE748701 


Homo sapiens CD48 molecule (CD48) 


9.68 


120 d 


A 

r\_ 


7n 


Dr\AAor\ 1 
_r WtoU I 


EE862627 


Homo sapiens chemokine ligand 26-like protein (CCL26) 


7.26 


120 d 




70 




tE805356 


Homo sapiens sorting nexin 10 (SNX10) 


5.93 


120 d 


A. 


.70. 


.P0001 26 


EE807554 


Homo sapiens guanylate-binding protein 4 (GBP4) 


4.97 


120 d 


A. 


.70. 


.P0 16471 


EE829459 


Homo sapiens N-myc upstream-regulated 1 (NDRG1) 


4.40 


120 d 


A. 


.70. 


.P039346 


MM 001009789 


i-ir\mr\ cnDiona nrjncitarilanHin F rpir&Titr^r l^-P\ fPTl^FR^ 
/ IDI 1 ID jUUICI Ij Ltl UjLuU IQJ 1 HI 1 1 1 1 trUcLJlLII \\ r j \T 1 \3\ \\) 


432 


1 20 d 


A. 


.70. 


.P011531 


EE776172 


Homo sapiens CCAAT/enhancer-binding protein (C/EBP), beta (CEBPB) 


4.14 


120 d 


A. 


.70. 


.P045436 


EE827672 


Homo sapiens glutathione S-transferase theta 1 (GSTT1) 


4.13 


120 d 


A. 


.70. 


.P025036 


EE746884 


Homo sapiens regulator of G-protein signaling 1 (RGS1) 


11.22 


135 d 


A. 


.70. 


.P049266 


EE815004 


Bos taurus mRNA for T-cell receptor alpha chain, leader sequence, variable and 
joining region (TRA2) 


7.06 


135 d 


A. 


.70. 


.P030156 


NM_001 009735 


Homo sapiens ISG15 ubiquitin-like modifier (ISG15) 


4.28 


135 d 


A. 


.70. 


.P030132 


NM_001 104930 


Homo sapiens receptor (chemosensory) transporter protein 4 (RTP4) 


4.20 


135 d 


A. 


.70. 


.P001016 


FE027322 


Homo sapiens eukaryotic translation initiation factor 2, subunit 2 beta (EIF2S2) 


3.52 


135 d 



Differentially expressed probes were detected from five f-test analyses using GeneSpring 10.0 (Agilent Technologies). 

1 . Original probe name in the microarray analysis. 2. The GenBank Accession Number of expressed sequence tags provided by the manufacturer according to the 
specific probe name. 3. The highly homologous human sequences acquired by a BLASTN search in NCBI (search setting: query coverage not less than 50%, E- 
value less than 1.00E-100). 



via leptin in the hypothalamus [65,66]. LEPR expression 
was higher in U70 than in T70. Whether the higher 
LEPR expression at prenatal stages is associated with 
high postnatal adipose deposition remains to be 
investigated. 

Both hyperplasia and hypertrophy continued in skele- 
tal muscle at 85 d. FBP2, which encodes a gluconeogen- 
esis regulatory enzyme that catalyzes the hydrolysis of 
fructose 1,6-bisphosphate to fructose 6-phosphate and 



inorganic phosphate, was expressed higher in T85 than 
in U85. A previous study suggested that FBPase partici- 
pates in some nuclear processes during the development 
and regeneration of skeletal muscle [67,68]. Additionally, 
PRG3 expression was also higher in T85. Previous stu- 
dies have shown that PRG3 plays a role in some early 
aspects of skeletal myogenesis [69-72]. 

The expression levels of PIK3R1, TSHR, and PON1 
were higher in U85 and U100 than that in T85 or T100. 
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Table 8 Validation of differential expression by 
quantitative PCR (qPCR) 

Ratio of gene expression between two sheep 
breeds 1 



Gene 


Method 


T70/ 


T85/ 


T100/ 


T120/ 


T135/ 






U70 


U85 


U100 


U120 


U135 


EEF1A2 


qPCR 


0.23 


1.36 


4.17 


0.57 


0.11 




microarray 


8.92 


16.21* 


1.78 


0.78 


0.98 


TXNRD / 


qPCR 


131.73 


382.54 


35.06 


2.3 


1.16 




microarray 


1.38 


1.14 


5.20* 


3.15 


0.94 


INSIG1 


qPCR 


0.82 


1.65 


0.32 


0.16 


0.10 




microarray 


1.10 


1.11 


0.34* 


0.88 


0.87 


LEF-1 


qPCR 


6.64 


21.65 


0.51 


0.78 


0.16 




microarray 


0.52 


0.59 


0.47* 


0.50 


0.57 


TM0D4 


qPCR 


81.5 


714.32 


323.03 


3.7 


0.02 




microarray 


2.80* 


3.52* 


1.54 


1.15 


0.80 


TPM2 


qPCR 


0.29 


0.87 


1.48 


0.17 


0.84 




microarray 


0.48* 


0.60 


1.54 


3.03 


1.10 


MY09A 


qPCR 


0.14 


0.15 


0.25 


0.31 


0.12 




microarray 


0.49* 


0.61 


1.02 


0.75 


1.20 


TRRAP 


qPCR 


0.06 


0.30 


0.85 


1.11 


0.42 




microarray 


0.51 


0.35 


1.79 


0.20* 


0.26 


TNC 


qPCR 


0.21 


0.54 


5.12 


6.23 


10.46 




microarray 


0.88 


0.69 


0.41* 


0.62 


1.78 


PIK3R1 


qPCR 


0.61 


0.44 


0.23 


1.06 


0.29 




microarray 


0.51 


0.36* 


1.03 


1.92 


1.1-1 


ACTB 


qPCR 


1.79 


0.16 


0.73 


0.35 


0.85 




microarray 


1.09 


0.55 


0.48* 


1.39 


0.47 


MYL6B 


qPCR 


0.35 


1.0 


3.44 


0.84 


0.85 




microarray 


0.42* 


1.14 


0.88 


1.16 


1 .04 


MYH1 


qPCR 


0.84 


1.36 


0.21 


0.57 


0.56 




microarray 


1.89 


5.60* 


1.56 


1.19 


1.09 


CA5Q1 


qPCR 


0.51 


0.76 


0.56 


0.07 


2.79 




microarray 


2.16 * 


2.31 


1.23 


0.82 


0.72 


MSTN 


qPCR 


0.81 


0.71 


0.54 


0.47 


0.6 1 




microarray 


0.69 


0.92 


0.62 


0.56 


0.81 


1. The data in the rows are the ratio of the 2 AacT values for the qPCR in a 



Texel sheep sample relative to that in the Ujumqin sheep sample at the same 
developmental stage. The microarray row provides gene expression ratios in 
Texel sheep samples divided by those in Ujumqin sheep samples at the same 
developmental stage. ^Differentially expressed genes in the microarray 
analysis with a statistically significant difference at p < 0.05. 



Previous studies have demonstrated that these genes are 
involved in fat metabolism and adipocyte development 
[73-76]. 

PDK4, a gene that was upregulated in T100, is a key 
regulatory enzyme involved in switching the energy 
source from glucose to fatty acids in response to physio- 
logical conditions. The pyruvate dehydrogenase complex 
occupies a central and strategic position in muscle inter- 
mediary metabolism and is primarily regulated by phos- 
phorylation/dephosphorylation [77]. In addition, PDK4 
is significantly associated with intramuscular fat and 



muscle water content [78], indicating that PDK4 is 
involved in meat quality. 

TXNRD1 is a variant of thioredoxin reductase, an 
important selenoprotein that maintains cellular redox 
balance and regulates several redox-dependent processes 
during apoptosis, cell proliferation, and differentiation 
[79-81]. LEF1 is a transcription factor involved in the 
regulation of myogenesis. Loss of Lefl -mediated repres- 
sion results in an increased number of cells expressing 
Pax-7 and Pax-3, suggesting that Wnt signaling via Lefl 
acts to regulate the number of premyogenic cells in 
somites [82]. In our study, LEF1 expression was higher 
in U100 than that in T100, which is consistent with the 
rapid proliferation of myofibers in U100. Although tran- 
scription factor activity is predominantly manifested at 
the protein level, the difference in LEF1 expression level 
indicated a potential difference of LEF1 between breeds. 

OVAR-DQB1 (sheep MHC class II) was expressed 
much higher in T100 and T120, particularly in T120 
(Table 7), compared with U 100 and U120. Results from 
Karpati et al. indicated that the class I molecule may be 
involved in the fusion of myogenic cells during muscle 
regeneration [83]. Honda and Rostami demonstrated 
that the expression of class I antigens on muscle cells is 
not only immunologically modulated but also develop- 
mentally regulated, and that these antigens may play a 
role in cell recognition and interaction during the myo- 
genic fusion process. The presence of the antigens, how- 
ever, was transitory, and they disappeared as myoblasts 
fused and differentiated into multinucleate myotubes 
[84]. 

CEBPB, an important transcription factor, acts as a 
indispensable regulator of adipocyte differentiation dur- 
ing adipogenesis [85-87]. We found that CEBPB expres- 
sion was higher in skeletal muscle in T120 than in 
U120, indicating a key difference in adipogenesis 
between Texel and Ujumqin sheep. However, which 
breeds contain more adipose mass than the other at this 
developmental stage remains to be investigated. 

TRRAP, which was downregulated in T120, in contrast 
to that in U120, is vital for embryonic survival and con- 
trol of the mitotic checkpoint [88-91]. Deletion of 
TRRAP leads to a reduced level of beta-catenin ubiquiti- 
nation, a lower degradation rate, and accumulation of 
the beta-catenin protein, whereas TRRAP knockdown 
results in abnormal retention of beta-catenin at chroma- 
tin and concomitant hyperactivation of the canonical 
Wnt pathway [92]. However, the canonical Wnt path- 
way is downregulated in the absence of myostatin 
through beta-catenin, whereas the Wnt/calcium pathway 
is upregulated [57], suggesting that TRRAP may nega- 
tively regulate skeletal muscle development between the 
canonical Wnt and the Wnt/calcium pathways through 
beta-catenin at 120 d. Further efforts are warranted to 
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test the effect of TRRAP on skeletal muscle 
development. 

Other than the genes discussed above, we found that 
MSTN expression was lower in Texel than in Ujumqin 
sheep through the five development stages (Table 8). 
This result indicated that MSTN downregulation con- 
tributed to the hyperplasia and hypertrophy of myofibers 
in Texel sheep through the second half of gestation. 
Moreover, the genes discussed above are potential 
myostatin targets for further investigation. 

Conclusions 

We demonstrated that fetal myofiber number increased 
sharply in a pulsed-wave mode at a particular develop- 
mental stage and was not augmented evenly across 
developmental stages in large mammal. The two surges 
in myofiber hyperplasia occurred at 85 and 120 d in 
Texel sheep, whereas a single surge appeared at 100 d 
in Ujumqin sheep during the second half of gestation. A 
myostatin mutation changed the gene expression profile 
in prenatal skeletal muscle, particularly disrupting some 
pivotal signaling pathways governing muscle develop- 
ment and function at some developmental stages, which 
explains much of the variation in myofiber phenotypes 
between Texel and Ujumqin sheep. Further studies on 
the crucial DE genes and signaling pathways involved 
would be helpful for revealing the mystery of muscle 
development in mammals. 

Methods 

Animals 

All experimental and surgical procedures were approved 
by the Biological Studies Animal Care and Use Commit- 
tee, Shanxi Province, Peoples Republic of China. 
Seventy-eight Ujumqin and 54 Texel ewes were pre- 
pared. These purebred female animals involved were 
selected based on their age (3-5 years old), body weight 
(50-55 kg), and body size. After these animals were sub- 
jected to pre-feeding for 45 days, estrus was synchro- 
nized in all ewes using an implanted controlled internal 
drug release device (Pharmacia & Upjohn Pty Limited, 
Parramatta City, NSW, Australia) and intramuscular 
injections of pregnant mare serum gonadotropin 
(Ningbo Renjian Pharmaceutical Co., Ltd, Zhejiang, 
China), according to the manufacturers' protocol, fol- 
lowed by artificial insemination using the corresponding 
breed sire's sperm. Three pregnant ewes from each 
breed were subject to caesarean section to collect the 
fetuses at 70, 85, 100, 120, and 135 d of gestation, and 
then the 12 different anatomic skeletal muscles from 
each fetus were dissected and weighed (longissimus 
dorsi, semitendenosus, semimembranosus, gluteus med- 
ius, femoral quadriceps, gastrocnemius, serratus ventralis 
thoracis muscle, biceps femoris, adductor, supraspinatus, 



infraspinatus, and triceps brachii). Two samples were 
taken from each muscle at standardized anatomical 
sites. The first sample (3 g) was dissected and rapidly 
frozen whole in isopentane chilled over liquid nitrogen 
for histological examination. The second sample (5 g) 
was snap-frozen in liquid nitrogen for gene expression 
analysis. 

Histology analysis 

To examine the development of skeletal muscle in both 
sheep breeds, we measured the number and the dia- 
meter of the longissimus dorsi muscle fiber in Ujumqin 
and Texel sheep during five development stages using 
hematoxylin and eosin staining. Serial cross sections of 
10 urn thickness were cut at -20°C using a cryostat. 
Photomicrographs of five to eight randomly selected 
areas were used to estimate the the diameter and the 
total number of myofibers in the cross sections of mus- 
cle examined. Muscle fiber measurements were con- 
ducted using DT2000 (V2.0) image analysis software 
(Nanjing East-image Ltd., Nanjing, China). 

The primary IFS antibodies were rabbit-laminin anti- 
body (Sigma, St. Louis, MO, USA) and mouse monoclo- 
nal anti-pax7 (Developmental Studies Hybridoma Bank, 
University of Iowa, Iowa City, IA, USA). The corre- 
sponding secondary antibodies were anti-mouse immu- 
noglobulin (ALEXA488) and anti-rabbit 
immunoglobulin (ALEXA594) (Invitrogen, Carlsbad, 
CA, USA). Staining was performed according to the 
manufacturer's instructions. The IFS slides were ana- 
lyzed using a system of epi-fluorescence microscopes 
(TE2000U; Nikon, Tokyo, Japan). The satellite cell (SC) 
nuclei were stained green and the basal lamina, red. 
Images were subsequently analyzed using Image Pro 
Plus (v5.1) image analysis software (Media Cybernetics, 
Inc., Bethesda, MD, USA). The numbers of SCs and 
myofibers in each section were counted, and the SC 
concentration within a given area of view was also cal- 
culated for each slide. 

Data between the two breeds were analyzed using the 
t-test procedure in SPSS 13.0 (SPSS, Inc., Chicago, IL, 
USA). A p-value of < 0.05 was considered statistically 
significant. 

Microarray experiment 

We made five independent comparisons at the same 
development stages in both breeds using a i-test and 
detected 207, 115, 90, 151, and 72 DE probes with a 
greater than twofold change at 70, 85, 100, 120, and 135 
d of gestation, respectively. The new specialized and 
standardized transcriptome-wide sheep microarray (Agi- 
lent Technologies) contained 15,008 sheep probes 
sourced from RefSeq and UniGene. The oligo microar- 
ray was 15 K in size. The longissimus dorsi muscles 
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from ovine fetuses were investigated for gene expres- 
sion. The microarray of each development stage for 
each breed contained three biological replicates (n = 3) 
except the T100 team (n = 4). Thus, 31 hybridizations 
in total were performed. Data were deposited in the 
NCBI database (GSE23563). 

Total RNA was extracted from the frozen longissimus 
dorsi muscles using an RNeasy Mini Kit with DNase 
treatment (Qiagen, Valencia, CA, USA) according to the 
manufacturer's protocol. The total RNA from each indi- 
vidual was labeled with Cy3 fluorescence. The generated 
targets were mixed and subjected to hybridization with 
the Sheep Gene Expression Microarray (Agilent Tech- 
nologies), according to the manufacturer's protocol. 

Scanning of the microarrays was performed using a 
DNA microarray scanner (Agilent Technologies). The 
scanning parameters were a 5 um scan resolution, PMT 
with 100%, and 10% once, respectively. All probe sets 
whose signals were reported as "absent" during the flag 
call in all 31 arrays were filtered out, and those from 
probe sets called "present" on at least one chip were 
used for further analysis. Scanner output image files 
were normalized and filtered using Feature Extraction 
Software vl0.5 (Agilent Technologies). 

Statistical analysis of the microarray data 

The raw signal intensity data were normalized using the 
quantile method, and the data were transformed to a 
base-2 logarithm for further statistical analysis. To 
explore the difference in gene expression between 
breeds at each developmental stage, we conducted a 
comparative analysis with the unpaired unequal variance 
t-test (Welch's £-test) method in GeneSpringlO.O (Agi- 
lent Technologies). Five t-test contrasts were performed 
independently, based on the premise that no relation- 
ships existed between the data at different time points. 
The thresholds for significance were p < 0.05 and a fold 
change of twofold or more with FDR < 0.4. 

Gene ontology and functional annotation analysis 

Only 1508 of the 15,008 probe set were definitive in the 
Agilent Sheep Gene Expression Microarray. To further 
clarify the function of the DE genes in this study, we 
conducted a BLASTN search for highly homologous 
human sequences (search setting: query coverage not 
less than 50%, E value less than 1.00E-100). DAVID 
Bioinformatics Resources (http://david.abcc.ncifcrf.gov) 
[93,94] and IPA (http://www.ingenuity.com) [95] were 
both used for further functional analysis to identify 
which molecular events or cascades were involved dur- 
ing skeletal muscle development. Significance was 
expressed as a p-value, which was calculated using the 
EASE score (p-values < 0.05 were considered 
significant). 



Quantitative real-time PCR analysis 

Total RNA was used to make cDNA using the Prime- 
Script® 1st Strand cDNA Synthesis Kit (Takara Bio, 
Shiga, Japan). q-Real-time PCR was performed on an 
ABI 7500 instrument (Applied Biosystems, Foster City, 
CA, USA) using the Fast EvaGreen® Master Mix (Bio- 
tium, Hayward, CA, USA). Thermal cycling consisted of 
an initial step at 95°C for 10 min followed by 40 cycles 
at 95°C for 30 s and 62°C for 30 s. Fifteen DE genes 
identified in the microarray experiment were selected 
for validation by q-real-time PCR; RpLPO was chosen as 
the reference gene for normalization of all data because 
it was expressed more stably [26-29]. Each q-real-time 
PCR reaction (in 20 ul) contained 10 ul of Fast Eva- 
Green® Master Mix, 0.5 ul of each primer, 2 ul of nor- 
malized template cDNA, and 7 ul water. The q-real- 
time PCR measurements were performed in triplicate 
on each cDNA sample (n = 31), and gene expression 
was quantified relative to RpLPO expression using the 2- 
AACt method. Gene expression ratios were normalized 
to the RpLPO gene in the same sample. 

Additional material 



Additional file 1: Quality control and validity of the gene 
expression microarray. This file contains the analytic settings using 
GeneSpringlO.O (Agilent) and the detectable rate of Agilent Sheep Gene 
Expression Microarray. 

Additional file 2: Differentially expressed (DE) genes in skeletal 
muscle between Texel and Ujumqin sheep at 70 d. For Additional 
files 2-6, the f-test was employed to identify the DE genes between 
Texel and Ujumqin sheep at each developmental stage using 
GeneSpring 10.0 (Agilent). Only those genes with p < 0.05(p-value in 
Colum B) and an absolute twofold change or more(absolute FoldChange 
in Colum C) were considered DE genes for further analysis. In the 
"regulation" column, "up" or "down" indicates that the genes were 
upregulated or downregulated in Texel fetal sheep, respectively, 
compared with those in Ujumqin fetal sheep. 

Additional file 3: Differentially expressed (DE) genes in skeletal 
muscle between Texel and Ujumqin sheep at 85 d 

Additional file 4: Differentially expressed (DE) genes in skeletal 
muscle between Texel and Ujumqin sheep at 100 d 

Additional file 5: Differentially expressed (DE) genes in skeletal 
muscle between Texel and Ujumqin sheep at 120 d 

Additional file 6: Differentially expressed (DE) genes in skeletal 
muscle between Texel and Ujumqin sheep at 135 d 
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